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ABSTRACT 



We present a survey of nonradial mode properties in evolutionary sequences of RR 
Lyrae star models. Attention is focused on the modes that may be driven by the opacity 
mechanism and on those that may be excited as a consequence of the 1:1 resonance with 
^-j- ■ the radial pulsation. 

Qualitatively, all the models share the same properties of the nonradial modes. At the 
quantitative level, the properties are to a large extent determined by the radial mode pe- 
riods. There is only weak dependence on the star metallicity and no apparent dependence 
on the evolutionary status, that is on the helium exhaustion in the convective core. 
On . In the whole range of RRab and RRc star parameters we find unstable nonradial 

modes driven by the opacity mechanism. An instability of radial pulsation to a resonant 
excitation of nonradial oscillations is also a common phenomenon in both types. We 
l' discuss a possible role of nonradial modes in amplitude modulation observed in certain 

RR Lyrae stars. 
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1 Introduction 



Whether or not nonradial modes play a role in RR Lyrae pulsation has been 
a matter of speculation for some time. Recently, Olech et al. (1999) pre- 
sented first circumstantial evidence for nonradial modes presence in three 
RRc variables in M55. The evidence was based on the power spectra which 
revealed presence of additional modes whose frequencies could not be at- 
tributed to radial modes. Similar power spectra were subsequently found in 
several other RRc anc RRab stars ( Olech et al. 1999b, Kovacs et al. 1999, 
Moskalik, P. 1999). 

Earlier (Kovacs, 1993; Van Hoolst and Waelkens, 1995) proposed the 
1:1 resonant excitation of nonradial modes in a radially pulsating star as 
an explanation of the Blazkho type modulation. Manifestation of this effect 
in periodograms is an occurrence of the equally-spaced side peaks around 
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the main frequency. Calculations provided by Van Hoolst et al (1998) con- 
firmed the plausibility of this idea. These authors studied stability of radial 
pulsation with use of the third amplitude equation formalism. They found 
that there is a high probability of a resonant excitation of a low ^-degree 
mode. However, their calculations were done only one stellar model. In 
this paper we apply the same formalism, with one additional simplification, 
to investigate stability of radial pulsation in a large sets of RR Lyrae star 
models. An outline of the formalism and the results are given in section 4.2. 

Still earlier (Dziembowski, 1977) instability of static models of RR Lyrae 
stars to certain nonradial modes has been demonstrated. The driving effect 
is the same as for radial pulsation. A linear instability, however, is not 
a sufficient condition for excitation. Nonlinear calculations are required 
to determine the ultimate outcome of the linear instability. Because of 
enormous numerical complexity, for the nonradial modes such calculations 
have never been done. The instability could be saturated with excitation of 
a single mode, which seems the most common situation among RR Lyrae 
stars. However, for instance among S Set stars typical is excitation of many 
modes. We do not understand why is it so. 

In section 4.1 of this paper we present a survey of the unstable nonradial 
modes in RR Lyrae stars and we discuss potential identifications of the 
modes detected by Olech et al. (1999). 

2 Evolutionary models 

All the stellar models adopted in the present investigation, have been com- 
puted adopting the latest version of the FRANEC evolutionary code, which 
includes several upgrades of the input physics. Major improvements are the 
opacity tables for the stellar interiors as given by Rogers <fe Iglesias (1992) 
and low-temperature molecular opacities for outer stellar layers by Alexan- 
der & Ferguson (1994). Both high- and low-temperature opacity tables have 
been computed by adopting the Grevesse (1991) solar chemical mixture. The 
equation of state is the OPAL one (Rogers et al. 1996), implemented in the 
temperature-density region not covered by OPAL, with the equation of state 
of Straniero (1988), plus a Saha Eos in the outer stellar layers (see Cassisi 
et al. 1998, 1999, for more details). As for the calibration of the superadia- 
batic envelope convection, the mixing length calibration provided by Salaris 
& Cassisi (1996) has been adopted. 

For the present work, we have computed Horizontal-Branch models for 
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Figure 1: Horizontal-branch evolutionary sequences. In SI and S2 sequences 
the ZAHB composition is characterized by Z = 0.001, Yjjb = 0.243 and that 
in S3 by Z = 0.0002, Y HB = 0.24. Stellar masses are 0.65, 0.67 and O.74M , 
respectively, in SI, S2, S3. The symbols mark the models selected for the 
pulsational analysis. 
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two different assumptions on the heavy element abundance: Z=0.0002 and 
0.001 namely. In both cases, an initial Helium abundance equal to Y=0.23 
has been adopted. All the HB models have as Red Giant Branch progenitor 
a structure with mass equal to O.8M . This means that when computing 
the Zero Age Horizontal Branch (ZAHB) models we have accounted for 
the evolutionary values for the size of the He core mass and the surface He 
abundance (Yhb) at the He ignition corresponding to a O.8M progenitor as 
provided by our own evolutionary computations for the previous H-burning 
phases. 

In Fig. 1, we show the selected evolutionary tracks in the H-R diagram. 
The symbols along each track indicate the models adopted for the following 
pulsational analysis. 

3 Linear nonadiabatic calculations 

Oscillation properties of the selected models were studied with the method 
developed by one of us (Dziembowski, 1977). Its recent updated descrip- 
tion may be found in Van Hoolst et al. (1998). For nonradial modes the 
equation of nonadiabatic oscillations are solved numerically in the envelope 
and matched to the asymptotic solution for g- modes, which is valid in the 
deep interior of RR Lyrae stars. The reason is that beneath the matching 
point the Brunt- Vaisala frequency is much larger than the oscillation fre- 
quencies. The Cowling approximation is assumed, which is well justified for 
the modes considered. The weakest point in the adopted method is the one 
related to the treatment of convective transport, whose Lagrangian pertur- 
bation is simply ignored. This is certainly a poor approximation but it is 
not essential for the main aim of this work because effects of convection on 
radial and nonradial modes are nearly the same.. 

As an introductory example we plot in Fig. 2 the growth rates, 7 = ^s(u>), 
and frequencies, / = ^.{uj)/2tt for modes at the selected degrees i for one of 
the models we chose for the pulsation analysis. The temporal dependence of 
oscillations is assumed in the form exp(— \uit). Effects of rotation has been 
ignored. Thus, each point represents 2^ + 1 normal modes. 

Let us note two types of unstable modes. There are isolated rapidly 
unstable modes with the growth rates 7 > 0.01 d _1 and sequences of modes 
with much lower 7's. In the former group we find radial modes and modes 
with 1= 6 and 10 which belong to the class of strongly trapped unsta- 
ble (STU) modes defined by Van Hoolst et al.(1998). These modes have 
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Figure 2: Growth rates and frequencies for modes of indicated spherical 
harmonic degrees in a selected model from S3 sequence. The model is char- 
acterized by the following parameters logT e g- = 3.822, log(L/L ) = 1.717, 

y c = o.i7 ' 
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no counterpart in the adiabatic approximation. In the interior the eigen- 
function of such modes are - to good approximation - described as inward 
propagating internal gravity waves with exponentially decreasing amplitude. 

The growth rate behavior in the sequences of low degree modes reflect 
the trapping properties of the acoustic cavity. Still at I = 1 even for the best 
trapped modes more than 80 percent of the kinetic energy is contributed by 
the g-mode propagation. The trapping effect is weaker at I = 2 and 3 but 
then it begins to increase. Note the sharp peak of 7 in the I = 6 sequence 
near the first overtone frequency. The STU modes occur always between 
the two best trapped ordinary modes. 

For the occurrence of STU mode a sufficient trapping in the evanescent 
zone separating p- and g-mode propagation zone is needed. In our selected 
mode the STU fundamental modes appear at £ = 8. With increasing £ 
they tend to Kelvin (f or surface) modes. The instability continues well 
above £ = 100. We hesitate to give the upper limit because of the increasing 
uncertainty due to our crude treatment of convection. Near the first overtone 
the STU modes begin at £ = 5 and end at £ = 15. 

4 Survey of pulsational properties 

The H-R positions of stellar models selected for this survey were shown in 
figure 1. The models cover various stages of the central helium burning. This 
can be seen in figure 3, where central helium content is plotted as function 
of the effective temperature. In the same figure, the periods of the first two 
radial modes are plotted. Solid symbols are used to denote linearly unstable 
modes. Second overtone is also unstable in some of our models. However, 
because three is no observational evidence for second overtone excitation in 
RR Lyrae stars, in present survey we consider only modes in the vicinity of 
the first two radial modes. 

4.1 Opacity-driven modes 

The general property of all models considered is that the trapping effect is 
weak in the £ = 2 — 4 range. In the vicinity of an unstable fundamental 
radial there are always unstable £ = 1 modes. In some models, like the 
one used in Fig. 2, there is a frequency range where the £ = 2 modes are 
unstable as well but with much lower growth rates. Rapid instability occurs 
only for the STU modes, which began in most of the models at i = 8. The 
trapping pattern near first overtone is similar to that near the fundamental 
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Figure 3: Central helium content (Y c ) in selected models is plotted as func- 
tion of their effective temperature in the left panel. In the right two panels 
periods of the radial fundamental and first overtone modes are given for the 
selected models. The empty symbols are used if the mode is stable 
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Figure 4: The relative frequency separation between most trapped (subscript 
" T" ) nonradial modes and the fundamental radial mode (top) . The relative 
(radial to nonradial) growth rates (bottom). The I = 1 mode correspond to 
the local maximum of 7. The £ = 10 mode is one of STU modes. 



mode. Again the most unstable are modes of the I = 1 degree and then 
STU modes which begin at i = 5 or 6. The main difference is instability at 
all low degrees. 

In Figs. 4 and 5 we show, respectively for the fundamental mode and 
first overtone ranges, the frequency distances to corresponding radial modes 
and the relative growth rates for most unstable 1 = 1 modes and for the 
selected STU modes. The most unstable i = 1 modes as well as all the 
STU modes have always higher frequencies than the corresponding radial 
modes. The two plotted parameters vary in rather narrow ranges and their 
values are determined by the radial mode periods. The dependence on the 
abundance (Z) is most easily seen in the distances of the STU modes. The 
dependence on the evolutionary status (Y c ), for which one should consult 
Fig. 3, is not recognizable. 

The growth rates of STU modes are similar to those of the radial modes. 
If the instability is saturated by one of these relatively high degree modes 
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Table 1: Amplitude-modulated RRc stars in M55 (from Olech et al.1999) 



Star P p [d] 


{fp fs}/ fp 


A s / S p 


V9 0.316 


-0.028 


0.19 


V10 0.332 


+0.004 


0.44 


V9 0.316 


-0.098 


0.57 



the star would appear as a nonpulsating object. At I = 5 — 10 the can- 
cellation of the opposite sign contribution would reduce the disc-averaged 
light amplitude to at most few millimagnitude. There is no firm evidence for 
occurrence of nonpulsating stars in RR Lyrae stars in the RR Lyrae domain 
of the H-R. The linear theory does not yield us a hint why radial modes are 
so much preferred over the STU modes by stars. 

The secondary peaks in the three amplitude-modulated RRc stars dis- 
covered by Olech et al.(1999) cannot be explained in terms STU mode ex- 
citations. In Table 1 we provide data on the distances between the primary 
and secondary peaks and the relative V-amplitudes. 

The secondary peak amplitudes are still too large and in two cases the 
frequencies are lower than those of the main peaks. Interpretation in terms 
of the t = 1 modes is more plausible though not free of difficulties. Also in 
this case the secondary peaks position present certain problem. However, 
the problem is not so essential because always we find unstable I = 1 modes 
on both sides of the radial modes. Furthermore, we do not have arguments 
why radial modes should always have higher amplitudes. 



4.2 Resonant modes 

The criterion for the instability of radial pulsation to excitation of a resonant 
nonradial mode may be written in the following form 

f5R\ 2 Id 2 + k 2 

where 5R is the amplitude of radius variations in radial pulsation; D denotes 
the frequency distance between the radial and the nonradial mode; k is the 
damping rate of the nonradial mode in the limit cycle of the radial mode; and 
C is the coupling coefficient. The criterion is from Van Hoolst et al.(1998) 
only the notation is different. 

Evaluation of the quantities occurring in the r.h.s. of Eq.(l) requires, 
in principle, nonlinear calculations, which we have not done. However, for 
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Figure 5: Similar to Fig.4 but for the first overtone vicinity. Here I = 7 is 
the chosen STU mode. 
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Figure 6: Properties of the resonant modes (closest to the fundamental 
radial mode). From top to bottom, the relative moment of inertia, the rela- 
tive frequency separation between consecutive modes, the fractional damp- 
ing rate are plotted against radial mode period for the fundamental mode 
vicinity. The parameters plotted are important for the resonant excitation 
of the nonradial modes. Modes with degrees £ = 1, and 5 or 6 are most 
likely excited. 
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a crude evaluation of the probability that a radial mode of specified ampli- 
tude is unstable to parametric excitation of a nonradial mode, we need only 
linear mode characteristic provided in Figs. 6 and 7 and certain coupling 
coefficients. These coefficients were evaluated by Van Hoolst et al.(1998) for 
the stellar model they selected. Here we rely on a simple scaling of their 
numbers which we will explain below. 

Let Vt{A) denotes the excitation probability at the radial mode ampli- 
tude A = 5R/R. Then, if effect of rotation are ignored, we have 

r o if a a c 2 < k 2 

Ve{A) = { Min (l, i^JS) if A *C 2 > k 2 (2) 

where Au> = uJ^ n -i — ^e,n denotes the frequency distance between consec- 
utive g-modes of degree £. Note that Auj/2 is the maximum frequency dis- 
tance between the radial and the nearest nonradial mode and that V A 4 C 2 — k 2 
is the distance at the onset of the instability. 

The coupling coefficients, C, for various radial & nonradial mode pairs 
in the model of RR Lyrae star was explicitly calculated by Van Hoolst et al. 
(1998) . From their numbers we found an approximate relation 

Ck,e = bk~Y 1 , (3) 

with 

b = 27 and h = \72dT 1 , 

where k denote radial mode order (here k = and 1 for the fundamental and 
the first overtone, respectively), / denotes mode inertia evaluated assuming 
the same amplitude at the surface. Our additional simplification consists in 
adopting the same bk values for all our models. 

Final simplification, which we adopted after Van Hoolst et al. (1998), 
is the assumption that k = — 7 9 , where -j g is the damping rate due to dis- 
sipation in the g-mode propagation zone. This seems well justified because 
we consider the situation when the opacity driven instability is saturated 
by the radial mode and the resonant nonradial modes have almost the same 
properties in the outer layers. Consequently, there should be also an exact 
balance between the driving and damping also for the nonradial mode. Then 
-j g is all what remains. 

The joint probability of the instability is given by 



V(A) = 1 - J] [1 - P e (A)} 

e 



(4) 
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Figure 7: Similar to Fig.6 but for the first overtone. In this case the I = 1 
and 4 are most likely excited. 

For rotating stars we have to consider modes of different azimuthal num- 
bers and evaluate probabilities Ve >m . The effect increases probability of the 
resonant instability (Dziembowski et al. 1988). 

Very much like Van Hoolst et al. (1998) we find the maximum of proba- 
bility of the resonant excitation of an I = 1 mode, both for the fundamental 
and the overtone radial modes and then an £ = 5 or 6 mode for the funda- 
mental and an t = 4 for the overtone. This is why we selected these ^-values 
in Figs. 6 and 7. In addition, there are data for £ = 2 modes. Excitation 
of these modes is less likely than the £ = 1 modes because of higher inertia. 
The data on h/Io are important for evaluation of effects of rotation and 
magnetic field on radial pulsation. We will not discuss these effects here. 

In Fig. 8 we present results of calculation of the excitation probabilities 
for modes of selected degrees which yield the dominant contribution to the 
joint probability. We chose A = 0.075 for the fundamental mode and A = 
0.025 for the first overtone. These values correspond to the mean amplitudes 
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Figure 8: Probability of instability of radial pulsation in the fundamental 
and first overtone mode relative to excitation of nonradial modes of indicated 
degrees. 
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of radius variations in, respectively, RRab and RRc stars (e.g. Jones et al. 
1988,Capricorni et al. 1989, Liu & Janes, 1990, Jones et al. 1992). The 
lower probability of the first overtone instability is a direct consequence of 
the lower value of A. 

The probability of the resonant instability in most cases increases with 
pulsation period. The exception is the instability of the fundamental mode 
to higher degree modes. In this case damping in the g-mode propagation 
zone plays an important role. The value of k increases with Pq and t. The 
increase reduces the chances for the instability and ultimately prevents it 
(see Eq. 2). 

At typical amplitude of RRab the probability of excitation of t = 1 
mode is between 0.25 and 0.5. This is not so different from the incidence of 
Blazkho effect which is estimated to be between 20% and 30%. The joint 
probability of the instability is always higher than 0.5 and close to 1 in most 
cases. However, excitation of modes with £ > 2 may not lead to amplitude 
modulation. The incidence of Blazkho effect amongst RRc stars is lower. 
Kovacs et al.(1999) who analyzed data on a large sample of RRc stars from 
LMC found the effect in 1.4% of the objects. Our analysis suggest lower 
chances for the first overtone instability than the fundamental mode but 
not in such a disproportion. 

5 Conclusions and discussion 

Our survey shows that all RR Lyrae star models share all qualitative proper- 
ties of the nonradial modes. There is always a large number of unstable low 
degree modes with frequencies close to unstable radial modes. However, ow- 
ing to higher mode inertia, for most of nonradial modes the driving rates are 
much lower than those for radial modes. The exceptions are the strongly 
trapped (STU) modes which begin with I degrees 7 to 10 (depending on 
the model) at frequencies somewhat above the fundamental radial mode 
and with t = 5 or 6 with at frequencies somewhat above the radial mode 
overtones. These modes are characterized by growth rates similar to radial 
modes. However, we argued that these modes are not likely candidates for 
identification of oscillation detected in some RR Lyrae stars. More likely 
candidates are the 1 = 1 modes. Their driving rates are by nearly an order 
of magnitude lower than radial modes but it is well known thatthe growth 
rate is not necessarily a good predictor of the finite amplitude pulsation. 
We found also that parameters which determine the chances of the exci- 
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tation of nonradial radial modes through the 1:1 resonance do not vary much 
over the range of RR Lyrae stars parameters. According to our estimate the 
excitation has a high probability. In fact some nonradial modes should be 
excited in majority of the RRab pulsators and in a significant fraction of 
( 30%) of RRc pulsators. The actual number should be greater because we 
ignored effect of rotation. Our crude estimate, which we did not detail here, 
shows that the effest is significant already at equatorial velocities of few 
km/s. 

Why then the incidence of the anomalous behavior among RR Lyrae 
stars is relatively rare? We should stress, that a significant amplitude mod- 
ulation is not automatically implied by the nonradial mode excitation. If the 
nonlinear interaction between radial and nonradial modes leads to a steady 
pulsation with constant amplitude then the presence of the nonradial mode 
will not be easy detectable. A Blazkho-type amplitude modulation may 
arise then only if the nonradial mode is not symmetric about the rotation 
axis and it is of low degree. In this case the Blazkho period is equal to the 
rotation period. Another possibility is a periodic limit cycle in which the 
amplitudes of the two modes vary intrinsically with the period determined 
by the nonlinear interaction. 

The ulitmate answer regarding the presence of nonradial modes in RR 
Lyrae stars may be expected only from spectropy. A signature of such modes 
should be searched in the line-profile variations. Thus, high-resolution spec- 
troscopic observations of amplitude modulated RR Lyrare stars are encour- 
aged. 
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